Can Resonant Oscillations of the Earth Ionosphere Influence 
the Human Brain Biorhythm? 
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Abstract 

Within the frames of Alfven sweep maser theory the description of morphological features of geomagnetic 
pulsations in the ionosphere with frequencies (0.1-10 Hz) in the vicinity of Schumann resonance (7.83 Hz) is 
obtained. It is shown that the related regular spectral shapes of geomagnetic pulsations in the ionosphere 
determined by "viscosity" and "elasticity" of magneto-plasma medium that control the nonlinear relaxation of 
energy and deviation of Alfven wave energy around its equilibrium value. Due to the fact that the frequency 
bands of Alfven maser resonant structures practically coincide with the frequency band delta- and partially 
rJieia-rhythms of human brain, the problem of degree of possible impact of electromagnetic "pearl" type resonant 
structures (0.1-5 Hz) onto the brain bio-rhythms stability is discussed. 

Keywords: Ionospheric Alfven resonator (IAR); ELF waves; cosmic rays; magnetobiological effect (MBE) in 
cell; brain diseases statistics 
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1 Introduction 

Lately Nobelist L. Montagnier's group has published 
three articles deeply challenging the standard views 
about genetic code and providing strong support for 
the notion of water memory [1-3]. In a series of del- 
icate experiments [1, 2] they demonstrated the possi- 
bility of the emission of low-frequency electromagnetic 
Extremely Low Frequency (ELF) waves from bacte- 
rial DNA sequences and the apparent ability of these 
waves to organize nucleotides (the "building" mate- 
rial of DNA) into new bacterial DNA by mediation of 
structures within water [4]. 

Without going into details of physical justification 
of quantum- field interpretation of these results 1 , let us 
emphasize one significant experimental result of this 
group. This result is related to stable detection of ELF 
waves (<7 Hz) from bacterial DNA sequences. Obvi- 
ously, if the result is reproduced in similar experiments 
of another research groups, the unique importance of 
this fact is difficult to overestimate in understanding 
of living matter essence. 

First of all, it concerns not only research of drasti- 
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cally new fundamental properties of spatial-temporal 
structure of eukaryotic genome, but also refers to 
studying of equally important issues that are related to 
exogenous nature of <7 Hz ELF waves impact directly 
onto a human brain and its biorhythms. 

It is well known that the brain neurons constitute 
different types of networks that interact by means 
of electrical signals. Neuron networks configurations 
comprise electrical circuits of oscillatory type. Elec- 
trical oscillations with different frequencies correspond 
to different states of brain. These oscillations could be 
detected by brain electroencephalogram. 

Numerous investigations have shown that electri- 
cal oscillations of different frequencies dominate in a 
healthy human being brain at its different states [6, 7]. 
Transitions between brain activities happen not con- 
tinuously, but only in discrete steps, from one level 
to another. The rest state corresponds to the stead- 
iest alpha-rhythm, with frequencies laying within the 
frequency band from 8 to 13 Hz. freta-rhythm with 
boundary frequencies 14-35 Hz corresponds to brain 
work. The slowest oscillations at frequencies 0.5-4 Hz 
are typical for c?e/ta-rhythm which corresponds to deep 
sleep. At last, i/ieta-rhythm with frequencies from 4 
to 7 Hz dominates in the brain if the state of nuisance 
or danger appears. 
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At the same time it is known [8-12] that the weak 
magnetic fields impact on biological systems is a sub- 
ject of the biophysics section called magnetobiology. 
It studies the biological reactions and mechanisms of 
the weak fields action. Magnetobiology is a part of 
a general fundamental problem of the biological effi- 
ciency of the weak and ultraweak physicochemical fac- 
tors, which operate below the biological defense mech- 
anisms threshold, and so may be accumulated on a 
subcelluar level. 

It is necessary to note here that there is no accept- 
able physical understanding of the way the weak mag- 
netic fields cause the living systems reaction [9] so 
far, although it has been experimentally found that 
such fields may change the biochemical reactions rate 
sharply in a resonance- like way [9, 13]. The physical 
nature of this phenomenon is still unclear, and it forms 
one of the most important, if not a general, problem 
of magnetobiology which includes the co-called "kT 
problem" . 

The problem consists in the fact that the weak mag- 
netic field energy (say, geomagnetic filed) of the same 
order as the kT heat energy is distributed over the 
volume 12 orders of magnitude larger (which approx- 
imately corresponds to a cell size). In such form the 
problem of the biological impact of the low-frequency 
magnetic oscillations has two aspects [9]: 

• what is a mechanism of the weak low-frequency 
magnetic signal transformation that causes the 
changes on the biochemical processes level of kT 
order? 

• what is a mechanism of such stability, i.e. how 
do such small impacts not get lost on the heat 
disturbances of kT order background? 

Not going into details of this complex and funda- 
mental problem, the essence of which is expounded in 
review by [9], let us note that in spite of the stated 
magnetobiology difficulties, there are serious reasons 
to believe that the main features of the magnetobio- 
logical effect are reliably established in numerous ex- 
periments and tests and are reproducible on differ- 
ent experimental models and under different magnetic 
conditions. On the other hand, the answers to the 
above-mentioned questions " lie" in the nonequilibrium 
thermodynamics field. "It is generally known that 
metabolism in living systems is a combination of pri- 
mary non-equilibrium processes. The origin and break- 
down of biophysical structures at time smaller than the 
time of thermalization of all degrees of freedom in these 
structures provide a good example of systems that are 
far from equilibrium where even weak field quanta can 
be manifested in system's breakdown parameters. In 



other words, if the life (thermalization) time of cer- 
tain degrees of freedom interacting with field quanta 
is larger than the system's characteristic time of life, 
then such degrees of freedom exist in the absence of 
temperature proper. Therefore, a comparison of their 
energy changes due to field quanta absorption with kT 
has no sense" [14]. The candidates for the solution of 
this problem today are the mechanism of the molecule 
quantum states interference for the idealized protein 
cavity [9, 14] and the mechanism of the molecular gy- 
roscope interference [9, 15]. 

Turning back to experiment, let us examine the pos- 
sible physical causes of the low-frequency geomagnetic 
fields generation and the consequences of their impact 
on the eucaryotic cells in short. 

It is well known that Earth's atmosphere between 
dense ionized shell called ionosphere (at an altitude 
of 100 km) and Earth's surface, possesses the elec- 
tromagnetic resonant properties (Fig. 1 [16]). Hence, 
resonances of the spherical cavity "Earth's surface 
- ionosphere" manifest themselves in electromagnetic 
quasi-monochromatic signals that permanently present 
nearby Earth's surface and has certain impact onto 
the Earth's biosphere. Among resonances of this type 
in the frequency band between (0.1-10) Hz the most 
known and studied is the so called Schumann reso- 
nance at the frequency of 7.83 Hz. This resonance 
is observed for electromagnetic waves with the wave- 
length exactly equal to the Earth's circle. Schumann 
resonance has drawn attention of physicians practi- 
cally immediately after its discovery in connection with 
studying of impact of electromagnetic radiation onto 
the a/p/m-rhythms of human brain which lie within the 
8-13 Hz frequencies band. 

Thunderstorms feed the Schumann resonator eter- 
nally. Initial frequency spectrum of electrical dis- 
charges (lightnings) during thunderstorms represents 
practically white noise. The resonant systems of near- 
Earth space filter out corresponding parts of the spec- 
trum which is shown in Fig. 2 [17]. 

Ionospheric Alfven Resonator (IAR) is being consid- 
ered along with Schumann resonator as a near-Earth 
resonant system, as well. In particular, with the help of 
IAR it was possible to explain new resonant radiation 
in the frequency band of 0.1-10 Hz [18]. This radia- 
tion was discovered in 1985 and characterized by quasi- 
periodic modulation (within frequency range of 0.5-3 
Hz) of the oscillations. This modulation appears above 
the background noise electromagnetic spectrum of the 
atmosphere and has regular daily variation (Fig. 2). It 
is not difficult to show [18], that within the frame of 
IAR model the resonant frequency f res of these oscilla- 
tions is defined by ionosphere layer thickness I, Earth's 
magnetic field strength H Earth > and concentration n of 
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Figure 1: Earth and surrounding electromagnetic resonant ob- 
jects (adapted from [16, 17]). a) Air gap at the altitudes of 0- 
100 km is the global Schumann spherical resonator with 7.83 Hz 
resonant frequency; the altitude region of 100-1000 km is a dense 
ionized shell (ionosphere). Inside its mass ionosphere Alfven res- 
onator with the first resonant frequency that varies in time within 
the limits of 0.5-3.0 Hz is located. Geomagnetic field lines lie 
above the ionosphere and are shown in red. High-energy protons 
cross this tubes. This Geomagnetic field line rests upon magneto- 
conjugated regions of the Earth's ionosphere which all together 
form the resonator of, so called, magnetosphere Alfven maser 
which generates "pearl" type electromagnetic signals. Traffic di- 
agram for the particles in radiation belt is also depicted in the 
figure. Particles with velocities being inside the loss cone (green 
line), possess small transversal velocities and fall into dense layers 
of the atmosphere, whereas particles outside the loss cone (blue 
line with arrow) possess bigger transversal velocities and are cap- 
tured by geomagnetic tube-trap due to their reflections from the 
magnetic mirrors of the ionosphere; b) Density of charged par- 
ticles in plasma of Earth's ionosphere versus altitude. 
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where va is Alfven velocity. According to [18] the spec- 
trum structure shown in Fig. 2 is defined by resonant 
frequency f res and its harmonics. For typical values of 
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H E arth ~ 0.4 E, M ~ 1.5 • 10— g, n 
I ~ 500 km this estimation gives f res ~ 2 Hz. Ap- 
parently, this estimation is in a good agreement with 
experimental frequencies estimations of the detected 
radiations 0.5-3 Hz [18]. 

Here it is interesting to mention another impor- 
tant role of IAR properties in affecting dynamics of 
larger scale resonator for Alfven waves: magneto- 
spheric resonator for Alfven waves Alfven Resonator 
(AR), formed by geomagnetic field line resting upon 
magneto-conjugated regions of Earth's surface. High- 
energy protons may cross geomagnetic field lines of 
the resonator and excite ultra-low frequency (ULF) 
electromagnetic oscillations practically in the same 
frequency band: 0.2-5 Hz , due to maser effect for 
the trapped protons and self-oscillatory mode of this 
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Figure 2: Electromagnetic noise spectrum structure for mid- 
dle latitudes has a pronounced resonant structure (adapted 
from [17]). In the daytime (c) and d)) the spectrum has a peak 
associated with Schuman resonance at 7.83 Hz, while at night 
time (a) and b)) electromagnetic noise produced by lightning's 
radiation at the frequencies below Schumann resonance is filtered 
out by Ionosphere Alfven Resonator (IAR) . In the years of solar 
activity maximum the noise spectrum at night time is similar to 
that of daytime. 



resonator [19]. This generator was called as mag- 
netospheric Alfven maser [16, 17, 20-22], which is 
schematically shown in Fig. 1. The signals generated 
via this mechanism are often referred to as "pearls". 
The spectral dynamical characteristics of the "pearls" 
and their temporal dependencies had been a mystery 
for researchers until recent times, a remarkable fact 
had been discovered recently consisting in the strong 
negative correlation between intensity level of low- 
frequency resonant lines in the atmospheric noise ra- 
diation spectrum and solar activity [17], and, corre- 
spondingly, between intensity of resonant radiations 
of "pearl" type and solar activity [23-25], that di- 
rectly correlate with predictions for IAR models [17] 
and Alfven maser ones [22]. 

Due to the fact that frequency bands of Alfven 
maser resonances practically coincide with the fre- 
quency band of delta- and partially i/ieia-rhythms of 
human brain, the problem of possible impact of electro- 
magnetic fields of "pearl" type onto stability of men- 
tioned brain bio-rhythms arises. 

Thus, investigation of possible direct correlation 
between the values of average annual frequencies 
of resonant electromagnetic signals of "pearl" type 
appearance, which have certain impact onto brain 
biorhythms, and rate of average annual mortality be- 
cause of diseases due to various abnormal functioning 
of human brain was the subject of this paper. Let us 
describe the basic laws of the electromagnetic fields 
generation and oscillations in the abovementioned res- 
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onators. 

2 Types of relaxation oscillations in 
Alfven maser 

Below we present a brief analysis of solutions for the 
equation describing small oscillations of the wave en- 
ergy around its equilibrium state in the Alfven sweep 
maser [20]: 

w + 2vw + 0|w = 0, (2) 

where 

E — Eq 
Eq 

and 

2v < n R . (4) 

Here Or and 2vr are characteristic frequency and 
oscillations decrement in Alfven resonator, respec- 
tively; w is the oscillation of the Alfven wave energy 
with respect to its equilibrium value Eq. The latter is 
defined for the case of no changes in the Earth's iono- 
sphere. These changes define reflection coefficients of 
Alfven waves and, consequently, their attenuation in 
the resonator under consideration; r rec is character- 
istic recombination time in ionosphere plasma; Nq is 
total number of fast particles in Geomagnetic field line 
having unit cross section at the ionosphere level; rtj s is 
electron concentration in the ionosphere; \ = {K^g), 
k is wave vector; v g is Alfven waves group velocity; 70 
corresponds to the steady state value of attenuation 
factor 

7 = |ln J R(w)|/r fl , (5) 

where R{lo) is a coefficient of Alfven wave reflection 
from the magnetic mirrors, that lasts over ionosphere 
and planet surface; r g is propagation time of elec- 
tromagnetic signal along geomagnetic field line be- 
tween magneto-conjugated regions of the ionosphere; 
$m = <Mw m )<5, 4>(uj m ) is normalized to unity Alfven 
wave amplification for one pass along the radiation belt 
(RB), while coefficient 5 equals to: 

4vre 2 /3 v 1 2 Q. pL 

= tt r , Po = — , W = -m e v , n A = -p, 

(6) 

where vq is typical velocity of particles, Q p l is the 
ion plasma frequency of background plasma in mag- 
netosphere RB equatorial cross-section; Q,^ is gy ro fre- 
quency in magnetosphere equatorial cross-section. It 
should be noted here, that the following approxima- 
tion for coefficient 5 m is used for taking into account 



the effect of the generated frequency sweep (drift): 

8 m = 5 + ^An is , (7) 

where 5q corresponds to equilibrium value of coefficient 
5, which is obtained for the case of stationary solution 
for differential equations which describe the dynamics 
of cyclotron instability in case of equilibrium electron 
density in the ionosphere. 

It is convenient for analysis of typical forms of re- 
laxation oscillations in Alfven sweep-maser to rewrite 
the equations (2) in the following way: 

w + Aw + TAw = 0, (8) 

with initial conditions 

w(0) = w , w(0) = 0, (9) 

where T = U 2 R /2v, A = 2v and w = [E{0) - E ] /Eq 
is a normalized initial energy of Alfven wave. 

The advantages of such representation of the Alfven 
waves energy relaxation oscillations become apparent 
during the study of the physical reasons for the time 
evolution of dispersion and, consequently, the morphol- 
ogy of such oscillations. For instance, it is easy to show 
that the Polyakov-Rappoport-Trakhtengerts equation 
(2) is equivalent to the following integro-differential 
equation: 

t 

w + (®^\ ■ 2v j e-^-^w^dt' = 0, w(0) = w . 


(10) 

As follows from (10), the medium memory function 
which characterizes its "elastic" properties, has the fol- 
lowing form: 

f(t-t')=u(t-t')-2ve- 2v( - t - t '\ t x = 1/2v, (11) 

where u(t — t') is the unit Heaviside function. Ob- 
viously, when v — > 00, equation (11) gains the 5- 
asymptotycs 

f(t-t')^S(t-t'), (12) 

while the equation (10) and, consequently, equation 
(2) as well, turns into a trivial relaxation equation of 
exponential type with the initial conditions: 

w = woe-( «/ 2 ") t = wo e - r ', (13) 

where 1/T = 2v/Q'j i is the time of w function "viscos- 
ity" relaxation to an equilibrium value wo, i.e. it is the 
relaxation time tm of the Maxwellian energy distribu- 
tion w to the equilibrium. 

Willing to preserve the properties of "viscosity" (i~m) 
and " quasi-elasiticity" (t\) of the medium in equation 
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(2) let us hereinafter consider the equation (2) in the 
form of (8) with any finite t\ = 1/A. 

So the characteristic equation, corresponding to (8) 
has the roots 



A: 



A 



1,2 



1 T \A - ^ 



r 

A 



rv A , (u) 



which are real when r\ < 1/4 so that the effective time 
of the medium aftereffect t\ < ta//4, where tm = 1/T 
is the time of the Maxwellian energy distribution set- 
tling. Particularly, in the case of a very short medium 
memory (rj -C 1) we have 



k\ = Ar?(l + t] + ... 
k 2 = A(l - V + . . . ) 



= r(i + 77), 

r(l-77)/r?. (15) 



In the case of ij > 1/4, which corresponds to Qr 3> 
2f^j, or Tjvf < 4ta (medium with a significant "elastic- 

ity") 

k\.2 = -A =F iw, w = -Aa/477 - 1. (16) 

Then the general solution satisfying the initial con- 
ditions has the form 



1 



w 



w 



w 



k\ + k 2 
E(Q) - E 



k 2 e~ klt 



k\e 



-k 2 t 



(17) 



In the case of i] < 1/4 the solution (17) takes on the 
following form: 



w 



w 



(1 + r/K 



-T(l+ V )t 



r/e 



-Xt 



V « 1, (18) 



which describes the exponential relaxation which is 
qualitatively different from e~ rt (a case of t\ = 0) only 
in the range < t < t\ = 1/A (Fig. 3a). Meanwhile 
for the case of 77 > 1/4 the solution (17) describes a 
new kind of mode 



w = woe 



-At/2 



sin(wi + if) 1 y/Ar] — 1 

: , r ? > 7> smy = 

sinp 4 

(19) 

in a form of a attenuated periodic relaxation (Fig. 3b). 

Then the expression (19), allowing for (3) and (17), 
may be represented in the following form 

E(t) = [25(0) - E ].e- xt / 2 Sin{ut + Lp) +E , E(0) > E, 

sin (p 

(20) 

where E(t) is the Alfven waves energy. 

It is known that the spectral analysis of experimental 
data corresponding to registration of magnetosphere 
radiation of "pearl" type or, in other words, geomag- 
netic pulsations Pel, allows one to reveal their inter- 
nal frequency structure [2], and the frequency inside 




Figure 3: Exponential (a) and oscillatory (b) types of relax- 
ations of Alfven wave energy perturbations 



of each "pearl" (separate packet of Alfven waves) in- 
creases from its beginning to the end [17]. Quanti- 
tative estimates of "pearl" parameters following from 
the theory above are in a good agreement with the re- 
lated experiments [17, 20]. Relying on that theory and 
experiment correspondence we will show below how 
the mentioned above morphological features find their 
explanations (within the framework of Alfven sweep- 
maser theory) on the basis of evolution of damping 
periodic oscillations (19) or (20) with taking into ac- 
count dispersion relaxation of Alfven wave energy to 
its equilibrium value. 

3 Maxwell distribution and relax- 
ation of Alfven wave energy dis- 
persion toward equilibrium value 

According to Eq. (17), deviation of Alfven wave energy 
fluctuation from average value is given by 



Aw(t) = w(t) - w(t) 



(21) 



where, according to (17) 

^) = W °^T^ {k 2 e- k ^ - k ie ~ k2t ) (22) 

and £(t) is normalized Gaussian noise with the follow- 
ing moments: 



m = o, mm = Ht-t') = <j)T5(t-t'). (23) 
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Let us consider the standard deviation 



2 

(Aw^ = J dt x J dhfl^-^-fae-^-^ 



t=\ 



_ kie -^ U )\ C{tl) {:(t 2 ). (24) 



Then, after integration (24) and taking into account 
(21) and (23) we obtain the following general expres- 
sion for dispersion of Alfven wave energy: 



(Aw) 5 



(k 2 -h) 2 



k 2 , 

ft 2 / i „-2fcii 



2fci 



2kik 2 
k\ + k 2 



1 



k\ , 



(25) 



Assuming that for t 3> tm = 1/T the energy distri- 
bution of Alfven wave is relaxing to Maxwell distribu- 
tion [26, 27]: 



(Aw)< 



t>i/r 



w 2 = e 2 



we have 

(l)T 



dw 
00' 



l (dw/d9) 



kT 

Eq ' 



(26) 



(k 2 - fci) 3 fc2/2fci - 2k 1 k 2 /(k 1 + fc 2 ) + k 2 /2k 2 ' 

(27) 

In the case of rj <C 1 (t\ -C tm) the dispersion re- 
laxation to its equilibrium value (26) is schematically 
shown in Fig. 4a. It is characterized by three relax- 
ation times: 



1 

2k 2 



+ k^V 
2 ki + k 2 



TA ' 2^ = 2T (1 - ??) - 
(28) 

In the case of rj > 1/4, when according to Eq. (16), 
k\ 2 = A/2 ± ia; and the relaxation character of dis- 
persion of Alfven wave energy to its equilibrium value 
(26) becomes an oscillatory one (Fig. 4b): 



(Aw) 



, dw 
~89 



-At 



1 - (l/V4r?cos(2wt + 3<p) 



1 - (l/v 7 ^?) cos3</? 



(29) 

where value is defined according to formula after 
Eq. (19) 

ip = arctan y/4rj — 1, (30) 

and the thermodynamical function dw/d9 by defini- 
tion is the thermal capacity Cy, which in the simplest 
case for plasma has the following form [26]: 



V 



{C ideal) v 



A 



2 ^3/2yl/2 ' 



const, 



(31) 



where C ideal is the thermal capacity of ideal gas. 




Figure 4: Relaxation of Alfven wave energy dispersion var(w) to 
equilibrium value for aperiodic (a) and oscillatory (b) character 
of relaxation in the medium with memory. 



Let us remind that the expression (29) taking into 
account (3) and (26) can be presented in the following 
form: 



(ae) 2 = a 



dE 



-\t 



1 - (l/^/4r/cos(2wt + 3v?) 



1 - (l/v 7 ^?) cos3y? 



(32) 

where 9j- = kT. 

Temporal behavior of average energy E (20) and 

^ l 2 

energy root-mean-square deviation (AE) (32) of 
relaxation oscillations of Alfven waves in magneto- 
plasma medium (medium with memory) are presented 
in Fig. 5. 

It is worth noting here, that such approach opens 
up nontrivial possibility for experimental numerical 
estimations of some important parameters of Alfven 
sweep-maser relaxation oscillations. For instance, the 
limit width of distribution (26) allows determining the 
plasma thermal capacity Cy for the given temperature 
9. In combination with the measured frequency of os- 
cillations £Ir it allows successively finding the values 
A(at t = it/uj), 7] (for any t < t\ = 1/A) and r m = 1/r. 

In other words, analysis of energy dispersion evolu- 
tion of Alfven sweep-maser relaxation oscillations al- 
lows finding experimentally (see Fig. 5) the values of 
plasma thermal capacity Cy, decrement A, relaxation 
time ("elasticity") t\ of magnetoplasma medium and 
settling time ("viscosity") tm of Maxwellian energy 
distribution (see Fig. 3). 



V.D. Rusov et al. o "Pearls" and Human Brain Biorhythm 



7 



'\> 

V i^t^ / -~- — i 




'V '/' »\\ 

'V '/' A* ' y?" - ^-— _ 

T '/' A v ' / v 




Ai 'A A"' / 2==- 
i\\ /// \ ^ — ' - 

A \ /// \ , ^ — 







particles with energy Wo ~ 200 /ce^) and allow for (6) 
and (35), we find that 



5s" 1 & f m ~ 1-25 ffz. 



(36) 



On the other hand, from (6), (35) and (36) it is not 
hard to find the value of Wo ■ 5, which (if we assume 
that rii S L ~ 3 • 10 2 cm -3 ) would be equal 



W Q 5 



10O| 



10 2 [. 



,-11 



(37) 



where the ionospheric plasma density nj S £ is defined 
Figure 5: Oscillatory relaxation of Alfven wave average energy E by the so-called plasma frequency 



and its dispersion var(E) in magnetoplasma medium (medium 
with memory). 



And finally, folowing [20], let us give some quanti- 
tative estimations. Accrording to [25], the recurrence 
period of the elements in "pearls" is about 50-300 s 
(Fig. 2), the bandwidth A/ fits into the range 0.05- 
0.3 Hz, and the dynamic spectrum tilt is 



dj_ 
dt 



2 • 10" 3 [Hz]. 



(33) 



In order to estimate the radiation parameters follow- 
ing from the sweep-maser theory [20] , let us consider a 
magnetic flux tube on the morning side of a magneto- 
sphere 2 at a distance of R « SREarth from the Earth 
center, where REarth is the Earth radius. In the frame- 
work of the sweep-maser theory it has been shown that 
the relaxation oscillations period, which characterizes 
the recurrence period of the elements in pulsations of 
the "pearl" type is 



Tr 



2vr 
Or 



2tt 



a-l 



Wn ■ 8 ■ 2So 



(34) 



where £Ir is the characteristic frequency of relaxation 
oscillations in Alfven resonator, a = Bm/Bi is the 
mirror ratio for the Earth's radiation belt, B m is the 
magnetic field at the ends of the magnetic trap, Bl is a 
magnetic field in equatorial section of magnetosphere, 
I is the effective length of a resonator, So is the equilib- 
rium precipitating protons flux density in the Earth's 
radiation belt. 

According to [20], the amplification curve 4>(co m ) 
(see. (3)) reaches its maximum when 



n 



1. 



(35) 



If we take into account the experimental values for 



10 2 



-1 



n A 



10 and A) ~ 2 • 10" 2 (for the 



2 Geomagnetic pulsations of "pearl" type are known [25] to 
appear primarily on the morning side of magnetosphere at mid 
latitudes at magnetically calm times. 



n 



P L 



y / Aire 2 n isL /m e , 



which determines the characteristic time scale of the 
plasma oscillations. 

Consequently, taking into account (34), (37) and the 
experimental data for a = 27 and I ~ R we derive 



Tr 



2tt 
Or 



2tt- 



10 4 



s, 



1/2 " 



(38) 



The period Tr obviously hits the experimentally ob- 
served range 50 -j- 300 s with the experimentally justi- 
fied value So ~ 10 5 cm~ 2 s~ l . 

Let us pass on to the dynamic spectrum tilt estima- 
tion: 

df_ _ df dn is ^ n df 
dt 



X ■ So 



(39) 



drii s dt drii s 

where rii s is the ionospheric plasma density, cm -3 . 

According to [20], the numerical calculation of 
df/drii s for the calm morning gives the value ~ 3 • 
10~ 6 cm 3 ■ s -1 . On the other hand there are reasons to 
believe that under weak magnetic storminess the pro- 



tons with energy about 200 keV (and x ~ 10 



cm 



- 1 ) 



flux density is So ~ 10 cm s . From this it follows 
that the dynamic spectrum tilt is 

df 



d[ 
dt 



X ■ So 



drii 



3 • 10" 



z ■ s 



(40) 



which corresponds to the experimental data [25] with 
a satisfiable accuracy. 

Therefore we may conclude that the Polyakov- 
Rappoport-Trakhtengerts sweep-maser theory [16-22] 
makes it possible to build a closed theory of genera- 
tion of a wide range of geomagnetic pulsations of the 
"pearl" type, which reside in the 0.1-^10 Hz band and 
are observed primarily at the mid latitudes under the 
conditions of a weak magnetic storminess on the morn- 
ing side of magnetosphere. In other words, it is shown 
that all the morphological features of geomagnetic pul- 
sations of the " pearl" type mentioned above find their 
adequate explanation in the framework of the Alfven 
sweep-maser theory. 
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Figure 6: Solar cycle (yellow) variations of geomagnetic pulsa- 
tions Pel (green) activity for about four (from 19 th through 21 st ) 
cycles [25]. 



4 On connections between varia- 
tions of geomagnetic pulsations, 
solar cycles and brain diseases 
mortality rate 

As mentioned above, the theory of formation of the 
"pearl" wave packets (geomagnetic pulsations Pel) in 
Alfven maser may help in solving another problem. 
This problem lies in the fact of strong inverse correla- 
tion between the appearance frequency of geomagnetic 
pulsations Pel and 11-year solar cycle. It has been 
found by means of long-term observations that geo- 
magnetic pulsations Pel activity is more intense (by 
factor of 10) during the periods of solar minima rather 
than in its maxima (Fig. 6). Below we will try to clar- 
ify this dependency within the frame of Alfven maser 
theory. 

Experimental observations of ionosphere have shown 
that steepness of electron concentration profile at the 
attitudes ~1000 km (see Fig. lb) is considerably de- 
creasing in the years of solar activity maxima [28] . This 
factor leads to decrease in the Alfven waves reflection 
coefficient from the upper layer of the IAR and, hence, 
to decrease in Q-factor of AR. Fig. 7 shows experimen- 
tal dependence of the reflection coefficient R of Alfven 
waves from IAR. It was obtained using ionosphere data 
for the minimum and maximum of solar activity [28]. 
One can see that appreciable decrease in the reflection 
coefficient R (and, therefore, worsening of the condi- 
tions for "pearl" generation in magnetospheric Alfven 
resonator (AR)) is observed for maximum of solar ac- 
tivity compared to its minimum. The explanation of 
such behavior of the reflection coefficient R and, con- 
sequently, behavior of the "pearl" generation rate is 
rather straightforward and is given below. 

The criterion for the wave generation in Alfven 




f, Hz 



Figure 7: Frequency dependence of Alfven wave reflection coeffi- 
cient R from ionosphere containing IAR [17]: 1 - for solar activity 
minimum, 2 - for solar activity maxim. Drop of plasma concen- 
tration at altitudes of 250-1000 km is much more pronounced in 
the solar activity minimum. It leads to the greater value of re- 
flection coefficient R (regions of R maximal values are indicated 
with arrows), where high rate generation of "pearls" takes place. 



maser according to (5) has the form [16, 17] 

R(u) -expTo > 1, (41) 

where Tq = ^r g is the logarithmic wave amplification 
at a single passing of AR (Fig. la), R{oj) is a coeffi- 
cient of Alfven wave reflection from the magnetic mir- 
rors, that lasts over ionosphere and planet surface. The 
" pearl" amplification changes little during the solar ac- 
tivity cycle and has its value considerably less than 
unity. Whereas the maximal value of reflection coeffi- 
cient R in the typical for the "pearls" frequency band 
of 0.2-5 Hz changes considerably according to Fig. 7 
and decreases in the years of solar activity maxima. 
So, it is follows from the (41) that the temporal vari- 
ations of IAR Q-factor affect the appearance rate of 
"pearls" generation. In other words, reflection coeffi- 
cient R behavior clearly explains (via the criterion of 
wave generation in Alfven maser (41)) the dynamics 
of "pearls" appearance rate, which, in turn, explains 
the reason for strong anticorrelation between "pearls" 
appearance and solar activity (Fig. 6). 

Using this dependence and having temporal evolu- 
tion of solar activity or, that the same, the tempo- 
ral evolution of Sun's magnetic field we may conclude 
about our principal knowledge of temporal evolution 
of "pearls" appearance over the past 100 years, at the 
least. Temporal dynamics of Sun's magnetic field is 
depicted in Fig. 8. Existence of strict anticorrelation 
between Sun's magnetic field and terrestrial magnetic 
field 3 (Y-component)[30] is seen from this figure, as 
well. 



Note that the strong (inverse) correlation between the tern- 
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Due to the fact that the frequency band of Alfven 
maser resonant structures practically coincides with 
the frequency band of deZta-rhythms and, partially, 
i/ieta-rhythms of human brain (see Fig. 7), the ques- 
tion naturally arises about the rate of possible influ- 
ence of global geomagnetic pulsations of "pearl" type 
onto stability of the above brain biorhythms. If such 
effect really exists, then one would expect positive cor- 
relation between variation of geomagnetic pulsations of 
"pearl" type Pel in the frequency band of 0.1-5 Hz and 
the death rate from disruption of brain diseases. Ob- 
viously, the choice in this case should concern only the 
currently incurable brain diseases so that their statis- 
tics were close to the real one and not being masked by 
intensive treatment. This fully applies to such diseases 
as malignant neoplasm of brain [31]. Their temporal 
dynamics in West Germany [31] is shown in Fig. 8. 
It is of interest that for our purposes the male statis- 
tics of infectious diseases (incl. Tuberculosis) in France 
[32] also applicable, that reflects, apparently, features 
of local spatial-temporal dynamics of magnetic field in 
Europe. 

It is easy to show, that degree of anti-correlation 
between temporal variations of Sun's magnetic field 
or, that the same, degree of direct correlation between 
frequency of "pearls" appearance and the number of 
deaths from considered diseases is high enough. This 
result is based on the experimental data on malignant 
neoplasm of brain [31], malignant brain tumor [36], 
brain lymphomas [37] and infectious diseases (incl. Tu- 
berculosis) [32]. In this way, according to Fig. 8, time 
arranged statistics of these diseases are lagging behind 
the variations of the Solar and Earth magnetic fields 
22-27 years and 10-15 years respectively. This delay 
effect on the one hand can be a consequence of the 
long-time hidden disease incubation period, but on the 
other hand opens up possibilities for prediction (at the 
time lag length) of behavior variations of indicated dis- 
eases by means of experimental observation of the ge- 
omagnetic field temporal variations. 

It is interesting to note here that a strong corre- 
lation between the galactic cosmic ray variations and 
cancer mortality birth cohorts has been discovered re- 



poral variations of magnetic flux in the tachocline zone and the 
Earth magnetic field (Y-component) are observed only for exper- 
imental data obtained at that observatories where the temporal 
variations of declination (dD/dt) or the closely associated east 
component (dY/dt) are directly proportional to the westward 
drift of magnetic features [29]. This condition is very important 
for understanding of physical nature of indicated above correla- 
tion, so far as it is known that just motions of the top layers 
of the Earth's core are responsible for most magnetic variations 
and, in particular, for the westward drift of magnetic features 
seen on the Earth's surface on the decade time scale. Europe 
and Australia are geographical places, where this condition is 
fulfilled (see Fig. 2 in [29]). 



cently [38, 39]. It was observed for population cohorts 
in five countries on the three continents. Previous ev- 
idence [39] has implicated a role for cosmic rays in US 
female cancer, involving a possible cross- generational 
foetal effect (grandmother effects). According to the 
assumption of the authors [38, 39], it may provide 
in-sight into the exploration of the role of germ cells 
as a possible target of this radiation and genetic or 
epigenetic sources of cancer predisposition that could 
be used to identify individuals carrying the radiation 
damage. And the conclusion about the galactic cosmic 
rays as a direct physical cause of the cancer mortality 
birth cohorts is based on a similar dependence of the 
total cancer age-standardized incidence rates and cos- 
mic ray rigidity from geomagnetic latitude (see Fig. 8 
in [38]). 

Not reducing the role of the physical mechanisms 
of radiation-induced effects formation and non-linear 
cell response in low doses of ionizing radiation (e.g. 
[40]), the study of which is a fundamental basis for the 
contemporary microdosimetry [41] , let us consider the 
possibility of indirect impact of electromagnetic res- 
onance structures (in 0.1-5 Hz band) of ionospheric 
Alfven maser on the cells of the birth cohorts through 
their direct impact on the germ cells of their parents. 
Fig. 8b, c shows a high level of inverse correlation be- 
tween the temporal variations of the solar magnetic 
field (or direct correlation between the "pearls" ap- 
pearance frequency) and cancer mortality rate of birth 
cohorts. 

Time lag between the inverse solar magnetic field 
and cancer mortality birth cohorts is ~6 years for UK- 
data and ~10 years for USA-data, as follows from 
Fig. 8b, c. At first sight it may seem to contradict the 
28- year lag between the galactic cosmic rays variations 
and cancer mortality birth cohorts, established in the 
paper by [38] basing on the data [42-44]. However, it 
may be explained by the known and hard-to-remove 
effect of the time shift in ice core data accompanying 
any 10 Be measurements (proxy of galactic cosmic rays) 
in ice cores of Greenland and Antarctica. On the other 
hand, theoretical verification of the actual 10 Be-data 
[45] and their comparison with the analogous data ob- 
tained from [42-44] indicates that the time lag between 
the galactic cosmic rays variations and cancer mortal- 
ity birth cohorts is ~6-10 years. 

There is also another more trivial justification of 
the 10-year time lag on the Fig. 8c. The variations 
of galactic cosmic rays are obviously a consequence of 
their modulation by the solar magnetic field. It means 
that magnetic fields of the solar wind deflect the pri- 
mary flux of charged cosmic particles, which leads to 
a reduction of cosmogenic nuclide (e.g. 10 Be and C) 
production in the Earth's atmosphere. In other words, 
cosmogenic nuclides (e.g. 10 Be and 14 C) are a kind 



V.D. Rusov et al. o "Pearls" and Human Brain Biorhythm 



10 




Figure 8: Time evolution (a) the variations of magnetic flux at the bottom (tachocline zone) of the Sun convective zone (see Fig. 7f 
in [33]), (b) fractional change in female breast cancer mortality for birth cohort in US (see Fig. 3b in [34]), (c) fractional change in 
female breast cancer mortality for birth cohort in UK (see Fig. 2b in [34]), (d) geomagnetic field secular variations (Y-component, 
nT/year) as observed at the Eskdalemuir observatory (England) [35], where the variations (dY/dt) are directly proportional to the 
westward drift of magnetic features, (e) Malignant brain tumor (brain stem) [36], (f) the number of deaths from ICD9 item n°191 
Malignant neoplasm of brain [31], (g) Brain lymphoma incidences in US [37] and (h) the mortality rates from infectious diseases 
(inch Tuberculosis) at ages 15-34 in France [32]. The curves (a) and (d) are smoothed by the sliding intervals of 5 and 11 years. 




1600 1700 1800 1900 2000 



Figure 9: Long-term cosmic rays reconstruction (from [45]). Cal- 
culated (grey curve) and actual annual 10 Be content in Greenland 
ice (dotted curve). Open circles represent the 8-year data from 
Antarctica [44]. Red line represents the 33- year moving average 
of the grey curve [45] . 

of a "shadow" of galactic cosmic rays on the Earth 
playing the role of a proxy for the solar magnetic vari- 
ability. Therefore the variations of the solar magnetic 
field and galactic cosmic rays (or 10 Be-proxy) must in- 
versely coincide which visually demonstrates the ex- 
perimentally justified result of the 1-year lag between 
10 Be and sunspot originally detected by Beer et al. 
[42]. After all, one could not expect anything else be- 
cause the galactic cosmic rays variations are caused by 
the solar magnetic field variations and not vice versa. 

Turning back to a direct physical cause of the cancer 
mortality birth cohorts, it should be noted that it is 
practically impossible to separate the possible radia- 



tive effect on germ cells (according to [34]) from the 
magnetobiological effect induced by such electromag- 
netic radiation as "pearls", since: 

a) electromagnetic radiation of the "pearls" type is 
generated as a result of the protons (a dominant 
component of the cosmic rays) passage through 
the Alfven maser resonator, which is a magneto- 
spheric magnetic flux tube resting upon the parts 
of ionosphere in the conjugate hemispheres of the 
Earth; 

b) the intensity variations of electromagnetic radia- 
tion of the "pearls" type - because of their origin 
- not only is correlated with the galactic cosmic 
rays variations, but also display a similar latitude 
dependence; 

c) magnetic field of the "pearls" can freely penetrate 
the human body just like any other magnetic field, 
because the human body tissues almost do not de- 
crease their intensity; indeed, the harmonic ampli- 
tude of the field with frequency w in a oscillatory 
circuit on the depth h inside the body is decreased 
f s times 

f s {u,h,a,n) = exp(-h/6), (42) 
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where the path till absorption 5 depends, according to 
[46], on the permeability fx (~ 1) and conductivity a, 
and is defined as follows: 5 = c{2-KfjLUJa)~ 1 / 2 . Since for 
uj < 10 6 we have S > 10 3 cm, for h ^ 10 cm from 
(42) we obtain the value f s ~ 1. 

Taking into account the stated properties and the 
known fact (e.g. [47, 48]) that the magnetic field may 
be a kind of an agent that amplifies the original cause 
(chemical impact or exposure to ionizing radiation) of 
the carcinogenesis, we may assume that the magne- 
tobiological effect induced by the electromagnetic ra- 
diation of the "pearls" type amplifies the cosmic rays 
radiative effect in the germ cells of the parents [34] . 

As one can easily see, the level of inverse correlation 
between the solar magnetic field variations (or direct 
correlation between the frequency of the "pearls" ap- 
pearance) and cancer mortality rate of birth cohorts 
is high enough. Also, according to Fig. 8c, the time 
series of this effect lag approximately 10 years behind 
the temporal variations of the Earth magnetic field. 
On the one hand, such delay effect may be a conse- 
quence of the cross-generational foetal effect [34], and 
on the other hand, it makes it possible to predict the 
discussed variations by experimental observation of the 
geomagnetic field variations. 

5 Conclusions 

In the frames of Alfven maser theory the description of 
morphological features of relaxation oscillations in the 
mode of geomagnetic pulsations of "pearl" type (Pel) 
in the ionosphere is obtained. These features are de- 
termined by "viscosity" and "elasticity" of magneto- 
plasma medium that control the nonlinear relaxation 
of energy and dispersion of Alfven wave energy to the 
equilibrium values. 

On the basis of analysis of the "pearls" generation 
criterion in Alfven maser (41) the physical reasons for 
strong anticorrelation of the appearance of "pearls" 
relatively to solar activity are discussed. Obviously, 
that a priori knowledge of the temporal evolution of so- 
lar activity or, that the same, the temporal evolution of 
Sun's magnetic field gives possibility to build the tem- 
poral evolution of "pearls" appearance over the past 
100 years, at least. The latter opens up a possibility 
for studying of positive correlation between "pearls" 
appearance rate and temporal variations of death rate 
from various disruptions of brain diseases. The anti- 
correlation rate between temporal variations of Sun's 
magnetic field or, that the same, direct correlation rate 
between of " pearls" appearance rate and the number of 
deaths from considered diseases was demonstrated to 
have the high enough value. This result is supported by 
the experimental data on malignant neoplasm of brain 



[31], malignant brain tumor [36], brain lymphomas [37] 
and infectious diseases (incl. Tuberculosis) [32]. 

The analysis of the known correlation between the 
galactic cosmic rays variations and cancer mortality 
birth cohorts observed for population cohorts in five 
countries on the three continents [34] let us suggest a 
hypothesis of a cooperative action of the cosmic rays 
and electromagnetic radiation of the "pearls" type on 
the germ cells of the parents which is responsible for 
the so-called cross-generational foetal effect [34] with 
a lag of ~6-10 years. 

In conclusion, we have obtained results clearly show- 
ing the possible impact of electromagnetic resonant ra- 
diations generated in ionospheric Alfven maser onto 
stability of human brain biorhythms, such as delta- 
rhythms and, partially, t/ieta-rhythms. 
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